Three-terminal thermoelectric transport through a molecular junction 
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The thermoelectric transport through a molecular bridge is discussed, with an emphasis on the 
effects of inelastic processes of the transport electrons caused by the coupling to the vibrational 
modes of the molecule. In particular it is found that when the molecule is strongly coupled to 
a thermal bath of its own, which may be at a temperature different from those of the electronic 
reservoirs, a heat current between the molecule and the electrons can be converted into an electric 
current. Expressions for the transport coefficients governing this conversion and similar ones are 
derived, and a possible scenario for increasing their magnitudes is outlined. 
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I. INTRODUCTION 

The investigation of thermoelectric phenomena in 
nanoscale devices at low temperatures has several in- 
teresting aspects. From the practical point of view, it 
is important to understand the heat flow and the dissi- 
pation because the heat generated by electric potentials 
used to switch-on transport currents inevitably induces 
decoherence in the quantum functioning of the device 
and also leads to dissipation. In bulk conductors, ther- 
moelectric transport necessitates an asymmetry between 
holes and electrons, which is usually small. In meso- 
scopic structures this asymmetry may be fairly high and 
can be also controlled experimentally. One would hence 
like to have a full picture of the symmetries and the inter- 
relations dominating the various transport coefficients of 
a small mesoscopic system, in particular the effects of 
inelastic processes. Indeed, when transport is through 
a molecular bridge, the tunneling electrons may undergo 
inelastic collisions with the vibrational modes even in the 
linear-response regime. This is because at finite tempera- 
tures the transport electrons may excite or de-excite the 
phonons and thus exchange energy with them. These 
inelastic processes modify the electronic transport coef- 
ficients, leading to the question of what, if any, are the 
analogues of the (bulk) Onsager-Casimir relations. An- 
other intriguing issue is the possibility to convert heat 
from the vibrations into an electric current between the 
electronic reservoirs, or vice- versa. 

Early studies of thermoelectric transport coeffi- 
cients of microstructures were based on the Landaucr 
approach^— which was also extended to include meso- 
scopic superconductors £ Once feasibility of measuring 
thermal and thermoelectric transport in atomic-scale 
samples had been established,— mesoscopic thermoelec- 
tric phenomena, e.g., peaks in the thermopower of a point 
contact orchestrated with the transitions among plateaux 



of the quantized conductance^ or oscillations (as a func- 
tion of a gate voltage) in the same coefficient measured 
on a quantum dot^— were detected and analyzed^ill 
The thermopower measured on nanotubes was found 
to be unexpectedly high, and this was attributed to 
a broken electron-hole symmetr y 12 ' 13 Similarly, nan- 
otubes exhibited enhanced thermal conductivity^ as did 
also silicon nanowiresJ^ The dependence of the ther- 
moelectric response on the length of the atomic chain 
has been recently computed within density-functional 
theoryJ£ Being based on the Landauer approach, the 
above-mentioned theoretical studies mainly focused on 
elastic processes of the transport electrons. Later on, 
effects of inelastic electron-electron processes and elec- 
tronic correlations (increasingly important at lower tem- 
peratures), as well as that of an applied magnetic field, on 
the thermopower produced in largeii and single-level^ 
quantum dots, and also in quantum wires^ were consid- 
ered. The effect of attractive electronic interactions on 
the thermopower was considered in Ref. |20| . 

Inelasticity of electronic processes should play a signif- 
icant role in thermoelectric transport through molecular 
bridges, also in the nonlinear rcgime^i Indeed, a density- 
functional computation of the nonlinear differential con- 
ductance of gold wires attributed changes in the I-V char- 
acteristics to phonon heating ; 22 ' 23 and the thermopower 
coefficient was proposed as a tool to monitor the excita- 
tion spectrum of a molecule forming the junction between 
two leads i 24 ' 25 It was suggested that the Seebeck effect 
in such bridges can be used for converting heat into elec- 
tric energy^ and to determine the location of the Fermi 
level of the transport electrons relative to the molecular 
levels, and also the sign of the dominant charge carriers, 
either for a molecular conductor ;2Zr— or for an atomic 
chain . 30 ' 31 This was confirmed experimentally: the See- 
beck coefficient as measured by STM yielded that in the 
case of the benzenedithiol family sandwiched between 



two gold electrodes the charge carriers are holes passing 
through the HOMO, whose location with respect to the 
metal Fermi level was determined from the magnitude of 
the coefficient £2- 

Inelastic electron-vibration interactions are not in- 
cluded in several of the theoretical studies devoted to 
molecular junctions (see, e.g., Refs. [Tfl [23, and [30) or 
are treated at off-resonance conditions, expanding them 
in the molecular-lead coupling^ When these interac- 
tions are ignored, the transport coefficients have the same 
functional form as in bulk conductors, with the energy- 
dependent transmission coefficient and its derivative re- 
placing the conductivity^ Notwithstanding the relative 
smallness, often, of the inelastic corrections to the ther- 
moelectric transport, their study is still of interest be- 
cause of fundamental questions related to the symmetries 
of the conventional transport coefficients, and since they 
give rise to additional coefficients connecting the heat 
transport in-between the electrons and the phonons. 




FIG. 1: (color online) A three-terminal system, modeled by 
a resonant level attached to two electronic reservoirs, hav- 
ing different chemical potentials and temperatures /il.r and 
Tl,r, respectively. An electron residing on the level interacts 
with its vibrational modes. The population of these phonons 
can be determined by the transport electrons (a "floating 
molecule") or by a coupling to a phonon source kept at tem- 
perature Tp. 

Here we study the heat and charge transport in a small 
mesoscopic (or nanometric) system depicted schemati- 
cally in Fig. Q] a molecule attached to two electronic (no 
phonons) reservoirs, held in general at different temper- 
atures, Tl,r, and at different chemical potentials, /iL,R- 
We distinguish between two (extreme) situations. In the 
first, the molecule is "floating" and is attached solely to 
the leads; then the vibration population is determined 
by the transport electrons alone. In that case, the sys- 
tem is a two-terminal junction. In the second case, the 
molecule is coupled to its own (typically a phonon-) heat 
bath which is kept at the temperature Tp, making the 
system a three-terminal one. It is then assumed implic- 
itly that the coupling of the molecule to that heat bath 
largely exceeds its coupling to the transport electrons. 
The latter is determined by our small parameter, the cou- 
pling between the molecule vibrations and the transport 
electrons, 7. Thus, we assume that the relaxation time 
due to the coupling to the heat bath, rp, is short on the 



scale r 2 /(7 2 wo) [see Eq. (JAI8I) in appendix [5], T being 
the level width on the molecule, due to the coupling with 
the leads and ujq is the frequency of the vibrations, h/rp 
may still be very small on all other physical scales, such 
as hioo and T. The phonon bath may be realized simply 
by an electronically insulating hard substrate (assuming 
that the large Kapitsa-type phonon thermal resistance 
between the lead and the sample is large enough to suffi- 
ciently reduce the thermal contact of the molecule to the 
substrate via the leads) , or a piece of such material touch- 
ing the junction, each of those held at a temperature Tp. 
A vacuum gap between the two separate substrates for 
the two leads would be ideal. However, with present fab- 
rication technology, this appears possible for a quantum 
dot but not for a small molecule. 

The consideration of the entropy production of such 
a three-terminal system is quite illuminating. Using the 
thermodynamic identity TdS — dE — [idN^ one finds 
that the dissipation at the left (right) reservoir leads to 
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(1) 



Here, — -El(r) is energy current emerging from the left 
(right) reservoir, while —Nub.) is t ne particle current 
leaving the left (right) reservoir. Adding to Eqs. ([T]) 
the entropy production of the phonon heat bath, Sp = 
Ep/Tp, where — Tp is the energy current leaving that 
bath, yields the total dissipation of the system, 

T P 

'-'P + ^L + "~>R = 7jT~ 

+ jr&- ^iA) + -£ (E R - Mr 7V r ) . (2) 
Charge conservation implies that 

N L + N R = 0, (3) 

while energy conservation requires 

T L + T R + T P = . (4) 

In the linear- response regime all three temperatures (see 
Fig. [J) are only slightly different, 

T L(R) - T ± — > 

T P = T + AT P , (5) 



and the chemical potentials differ by a small amount, 



Ml(r) =(i± 



Afj, 



(G) 



Expanding Eq. ([2]) and using Eqs. ([3]) and (@J yields 
S P + S L + S R =^f(-E P ) + ^I+^I Q , (7) 



where / is the net charge current flowing from the left 
reservoir to the right one, 



I=-^[^-N Rl 



(8) 



while Iq is the net heat current carried by the electrons, 
In = J w . - (u/e)I , with I w = --(& -En) . (9) 



'Q = j e ~ (p/e)I , with I E = - 2 ^l - ^R 



Finally, the heat current flowing from the phonon bath to 
the quantum system is simply given from the condition 
of energy conservation, 



-E P - E L + E R 



(10) 



Thus, the entropy production of our three-terminal sys- 
tem is a simple example of the general expressions for 
linear transport, consistent with the Onsager theory^ 

Since our molecular bridge is not necessarily at equi- 
librium within the transport process, it exchanges en- 
ergy with the phonons of the phonon reservoir by go- 
ing up/down in the vibrational ladder with absorb- 
ing/emitting a phonon in the bath. This is the physical 
origin of the current —Ep. On the other hand, when the 
molecule is floating, then — Ep vanishes. Since — Ep is 
proportional to the rate of change of the vibrational level 
population on the dot (see Appendix [5] for details) this 
in turn will determine the vibration population that will 
adjust itself according to the temperature and chemical 
potential differences applied to the electrons. In this sit- 
uation our device becomes a two-temrinal one, and the 
energy current carried by the electrons is conserved. 

In Sec. [IT] we outline our model, and give explicit ex- 
pressions for all three currents I, Iq, and —Ep, and in 
Sec. IHII we discuss them in the linear-response regime. In 
particular, we find there that by the three-terminal junc- 
tion, one may convert the heat current from the phonon 
bath into electric and heat currents carried by the elec- 
trons even at zero bias voltage and when Tl = Tr. In Sec. 
II VI we discuss the necessary conditions for this conversion 
to be established, i.e., the junction couplings to the elec- 
tronic reservoirs should not be spatially-symmetric and 
should not depend on the energy in an identical manner. 
We show that an opposite dependence on energy of the 
couplings to two electron reservoirs will tend to maximize 
the new transport coefficients we find. 



II. THE CURRENTS 

In our analysis, the molecular bridge is replaced by a 
single resonant level; when a transport electron resides on 
the level, it interacts (linearly) with the phonons. Such a 
model, which neglects effects of spin and electronic corre- 
lations, is applicable in the Coulomb blockade regime for 
low-energy molecular levels. (We also ignore the possi- 
bility of the Kondo effect to develop, namely the average 



temperature of the system should exceed the Kondo tem- 
perature.) The model Hamiltonian (see Fig. [1} is thus 



ri — 7lT ; + rLu + ri 



R • rl dot ~r '■-coup i 



rtf 



(11) 



in which Hl(r) is the Hamiltonian of the left (right) lead, 



^L(R) ~2^ 6 k(p) C k(p) C k(p) 

k( P ) 



(12) 



[using k(p) for the left (right) lead]. The Hamiltonian of 
the bridge, which includes the electron-phonon interac- 
tion, is 

^dot = e cta + w o( fotfo + 5) + 7(6 + fot )4c , (13) 

where luq is the frequency of the harmonic oscillator, and 
7 is its coupling to the transport electrons. (We use units 
in which h = 1.) Finally, the coupling between the dot 
and the leads is described by 

^coup = E(^ c Ic + he) + J2( V p4 c o + hc ) ■ (I 4 ) 



The operators c , c k , and c£ (c , c k , and c p ) create (de- 
stroy) an electron on the dot, on the left lead, and on the 
right lead, respectively, while tf (b) creates (destroys) 
an excitation of the harmonic oscillator, of frequency ojq. 
The electron distributions of the leads, /l and /r, are 
given by 

/l(r)M = (l + exP^LCR)^ - Ml(r))]) > (I 5 ) 

where /? L ( R) = l/(fc B T L(R )). 

The couplings of the leads to the resonance level broad- 
ens it. such that 



L(R) 



M 
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-HP)> 



(16) 



are the partial widths brought about by the left and the 
right leads. These couplings are treated to all orders, 
encompassing the case in which the transport electrons 
excite effectively the phonons (the dwell time of the elec- 
trons on the junction largely exceeds the response time 
of the oscillator, about Wg" 1 ), and also the inverse situ- 
ation. Strictly speaking, the Hamiltonian (fTT|) pertains 
to a "floating" molecule, which is not coupled to a heat 
bath of its own; however, the analysis presented above in 
Sec. Q] enables us to consider the three-terminal case (see 
Fig. [T|) as well. 

The explicit calculation of the currents is carried out 
up to second order in the electron-phonon coupling, us- 
ing the Keldysh technique^ and the details are given in 
Appendix [X] We find that the charge current consists of 
two terms, which can be related to clastic and inelastic 
transitions of the electrons through the junction [the first 



and the second terms in Eq. (fl7|) . respectively] 

1 = e J £|G5oMI 2 r L Hr R M(/ L H - / R ( 



to 



+ n-j 2^l%K)rlS5o(w-)r 

x (r R ( W+ )r L (o;_)F RL (u;) - L o r) , 
where we have introduced the abbreviations 

Wn 



VJ 



± 



Ul± 



-■() 



(17) 



(18) 



Here, Goo is the Green function of the dot, given by Eq. 
(|A5[) , and Qqq is its counterpart when the coupling to the 
phonons is ignored, i.e., 



ISSoMI 5 



i 



iT(u)/2 



(19) 



represents the bare Breit-Wigner resonance on the dot, 
with 



i» = r L H + r R ( w ) . 



(20) 



Finally, 



F aa ,(w) = N[l-f a (w + )]f „,((»_) 

-[l + #][l-/ a , («_)]/>+) (21) 

embodies the populations of the electrons (/l,r) and the 
phonons (N) . Note that the latter population is not nec- 
essarily given by the Bose-Einstein distribution; this is 
the case only when the molecule is strongly coupled to 
a heat bath of its own (this distribution is denoted be- 
low by N't). In the case of the floating molecule, the 
population N is determined by the transport electrons 
as explained below and in Appendix [X] . 

The energy current carried by the electrons, /g, [see 
Eq. (|9}] is shown in Appendix lAl to be 



h = 



^|GSoHI 2 ^r L Hr R H(/ L H - / R ( w )) 

|£|0So(« + )| a |65o(u;_)| a 

x (y [r R K)r R (u,_)*RRM - (R -► l)] 

+ w[T R ( W+ )r L ( W _)F HL (w) - (L o R)]) , (22) 



where again the first and second terms pertain to the 
elastic and inelastic contributions to the electronic energy 
current. The energy current carried by the phonons [see 
Eqs. (gj and flUJ)] is 



-Ef 



d'jj 
2^ 



%K)i a ia5o(«_)| s 



x 7 2 ^ Y, r Q ( w+ )r a ,( w _)F QQ ,H . (23) 

a, a'— L, R 

In the next section we examine these currents in the 
linear-response regime. 



III. THE LINEAR-RESPONSE REGIME 

The temperatures and the chemical potentials of the 
three-terminal junction are given by Eqs. ([5]) and ([5]). 
In the linear-response regime, one expands the currents 
[see Eqs. ©, ®, and (ID}] to first order in An, AT, 
and ATp . In order to express the resulting transport co- 
efficients in a convenient form, we note that all integrals 
resulting from the elastic processes include the function 



F e V) = /3/M[l-/H]|G 00 H| 5 



(24) 



where f(oj) is the thermal-equilibrium Fermi distribution 
of temperature T, and j3 = l/(fceT). The transport co- 
efficients coming from the inelastic processes include in 
their integral forms the function 

^ cl H=7 2 |^oK)| 2 |^o(^-)| 2 

xJV T /9/(w_)[l -/(«+)] , (25) 

where Nt is the thermal-equilibrium Bose distribution 
function of temperature T. 

The relations between the currents and the driving 
forces in the linear-response regime can be written in the 
matrix form 



(26) 



where the matrix of the transport coefficients, M., is 

G K X p 

K K 

X p 



I 




An/e 


J Q 


= M 


AT/T 


_-Ep. 




ATp/T 



M 



K p 



X* 



x p 
c p 



(27) 



Let us first describe the conventional transport coeffi- 
cients, pertaining to the transport by the electrons. In 
Eq. (|27| . G is the electrical conductance, 

G = G cl + G ino1 , (28) 

which consists of the contribution of elastic processes, 



G" 



e 

2^ 



d Wj F el ( w )r L ( w )r R (o;), 



(29) 



and the contribution of the inelastic ones 
p 2 

Qincl _ ' / J, , panel /■ 



=— / duF 1 
2tt 



>) 



r L( w +)r R ( w -) + r L( w -)r R ( w +) 



(30) 



Clearly, Eq. (|30l) corresponds to the two inelastic pro- 
cesses by which the transport electron excites or de- 
excites the phonon upon moving between the reservoirs. 
The transport coefficient yielding the thcrmopower and 
the Seebeck effect, K, and the one giving the main con- 
tribution to the electric thermal conductance, K2, also 
consist of two contributions each, 



K = K 
K 9 =K 



el 
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■1 



K 



Kind 
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(31) 



with 



K c 



2- 



du;F el ( U )(uj-n)T L (uj)T R (u;), 



T^inel _ / j, . ciinol 



= — j duF™ l (uj)(u; - „) 



x (r L ( w+ )r R ( w _) + r L ( w _)r R ( w+ )) , (32) 



and 



K5 



2vr 



d W F el ( w )( w - M ) 2 r L ( w )r R ( w ), 



K incl 



- — / dcuF^(Lu)(co - tf 



x (r L ( w+ )r R ( w _)+r L (a;_)r R (u; + )) 



(33) 



All other coefficients appearing in Eq. (|27|) result from 
the inelastic processes. One of them, KJj" , just augments 
the (conventional) ratio K2 between the heat current car- 
ried by the electrons and the temperature gradient AT 
across the junction, 



K p 2 = ^ I *****(«) 



x (r R ( w+ )r L (w_) + r L (« + )r R («_) 



(34) 



It therefore follows that the electron-phonon interaction 
just renormalizes slightly the conventional transport co- 
efficients of the two-terminal single dot junction, but does 
not lead to novel effects (see also Sec. IIVI below) . 

On the other hand, keeping the phonon bath to which 
the molecule is attached at a temperature different from 
those of the electron reservoirs leads to new thermoelec- 
tric effects. We find that there is an electric current flow- 
ing in response to the temperature difference ATp with 
the phonon bath, with the novel transport coefficient 



x P = e^o [ dwF in cl{w) 
2n J 

x (r R ( w+ )r L ( w _) - r L ( w+ )r R (u;_)) . (35) 



The same coefficient controls the heat current be- 
tween the junction and the phonon bath in response to 
the chemical potential difference between the electronic 
reservoirs. Likewise, there is a heat current flowing be- 
tween the electronic reservoirs in response to ATp , which 
is governed by a coefficient analogous to Eq. (|35l) . 



~p _I£q 

2tt 
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(36) 



r 



with the same coefficient governing the heat current from 
the phonon reservoir in response to the electronic temper- 
ature difference AT. Thus, the matrix of coefficients A4 
obeys the Onsager symmetry relations also in the three- 
terminal situation with the two types of carriers and their 
interaction. 

Finally, the coefficient C p gives the response of the 
heat current carried by the phonons to the temperature 
difference ATp, 



C l 



2tt 



dLuF mcl (io)T(uj + )r(uj_) 



(37) 



where we have used Eq. (f2T)|). 



IV. DISCUSSION 

Using a simple model, we have considered the ther- 
moelectric and thermal transport of electrons through 
a molecular bridge, in particular the subtle effects of 
the inelastic electron-vibrational mode processes. Of a 



paramount importance is the mechanism by which the 
vibration population is determined. 

When the molecule is not attached to any heat bath, 
the phonon population is determined by the voltage and 
the temperature difference across the junction. We show 
in Appendix lAl that in this case [see Eq. (|A16[) ] the heat 
current between the vibrations and the transport elec- 
trons is 



E, 



dN 
) ~dt 



(38) 



where N denotes the vibrational mode population. At 
steady-state that population does not vary with time, 
and consequently the heat current between the molecule 
and the junction vanishes. This requirement, in turn, 
fixes ATp in terms of A/i and AT, and consequently 
determines the vibration population [see Fig. [T]and Eqs. 
([5]) and ©]. In other words, the requirement that — Tp = 
yields 



X J 



,A/x 



X^ 



,AT 



= -C 



,AT D 



(39) 



and hence transforms the three-terminal junction into a 
two-terminal one, with 



I 

G - (X P J 2 /C P 
K - X P X P /C P 



K, 



K - X P X P /C P 

+ K P -(X P ) 2 /C P J 



Afi/e 

AT/T 



(40) 



In this situation we find that the inelastic processes mod- 
ify the transport coefficients, but do not give rise to any 
intriguing effects. 

On the other hand, when the molecule is attached 
(strongly) to its own thermal bath, see Fig. [TJ such 
that the system becomes a three-terminal junction, the 
vibrational modes and the transport electrons may ex- 
change heat, and a temperature difference between the 
phonons and the transport electrons can induce an elec- 
tron current between the electronic reservoirs. Likewise, 
a voltage between the latter can induce a heat current 
to the phonons. These two new transport coefficients, 
having two types of carriers and including inelastic pro- 
cesses, are related by Onsager symmetry. This situation 
is characterized by the appearance of new transport co- 
efficients that result solely from the inelastic transport 
processes [see Eqs. (f2l)|) and (|27p]. and requires the 
breaking of spatial symmetry between the two sides of 
the junction, IY ^ Tr. Note in particular the change 
of the relative sign of the combinations F R (w + )F L (w_) 
and r L (w + )r R (w_) between the expressions for the usual 
thermoelectric coefficients, Eqs. (j32l) . (|33| . and ([34]) , and 
the new three-terminal ones, Eqs. (|55|) and (|55)) . This 
change occurs because the latter expressions are for the 
heat currents from each lead to the phonons and not be- 
tween the two leads. The analysis of the above com- 
binations of the T's can tell us how to maximize the 
new, three-terminal, thermoelectric coefficients. Usually 
the w-dependence of the resonance widths T's is not too 
strong. Let us then expand them around the running w, 



r L (u/)=r L (u;) + ,4 L (u/-u/)- 



with an analogous expansion for T R . The crucial quantity 
is the one in parentheses on the right-hand side of Eq. 
(|55|) . To order uj it gives 



u;o(A R r L (u;) - A l T r (uj)) 



(41) 



To increase the usual thermopower, we want the trans- 
mission to depend strongly on energy. In our case, to 
make the two terms in Eq. (|4i"j) add and not tend to 
cancel, we also want T L and T R to have opposite de- 
pendencies on the frequency. One way to effect this is 
to have a lead with an electron-band material on the left 
lead, and one with a hole-band material on the right lead. 
This will however decrease the values of the usual two- 
terminal thermal and thermoelectric coefficients. Hence, 
more down-to-earth estimates of the new thermoelectric 



coefficients require realistic descriptions of the molecu- 
lar bridge, which will depend on the type of molecules 
involved and other parameters of the system. 

In order to elucidate the above considerations, we com- 
pute the coefficient governing the conversion of heat from 
the phonon bath into a voltage difference across the 
bridge, 



~.p _ 



eX p 
TG ' 



(42) 



where X p is given by Eq. (|3"5j) . (This definition follows 
the conventional one for the thermopower.) Let us as- 
sume that the left reservoir is represented by an electron 
band, such that the partial width it causes to the reso- 
nant level is given by 



r L («) 



L 



(43) 



u v - oj c 
while the right reservoir is modeled by a hole band, with 



r R (w) = f 



(44) 



Here, w c is the bottom of the conductance band (on the 
left side of the junction), while u> v is the ceiling of the 
hole band (on the right one). The energy integration 
determining the various transport coefficients is therefore 
limited to the region uj c < uj < uj v . (For convenience, we 
normalize the L's by the full band width, u v — uj c .) 
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FIG. 2: (color online) The coefficient S , Eq. (|42|>. as a func- 
tion of /3tv for /3r L = 0.2 (thin line), 0T h = 1 (dashed line), 
and /3Fl = 5 (dotted line). Upper panel, Fr = Tl, lower 
panel T R = 0.7r L [see Eqs. (J43J) and flU)]. The total band- 
width is determined by /3u; c = —fiujv — 100. 

Measuring all energies appearing in the explicit ex- 
pressions in units of the temperature /3 _1 , and choosing 



jU = eo=0 for simplicity, we obtained the curves shown in 
Figs. [5J One observes that the magnitude of the effect 
is non monotonic in the value of the vibration frequency 
loq: being the outcome of inelastic processes, it vanishes 
at ujq = 0, and also at /3wo 3> 1, since then the vibrational 
level population becomes very small. It increases for val- 
ues r L ( R ) which are smaller than the temperature, and 
also increases when IY ^ r R , because then the electric 
conductance G becomes smaller. 



Appendix A: Details of the currents' calculation 

The particle and the energy currents can be expressed 
in terms of the electronic Kcldysh Green functions, in 
particular the Green function Goo on the dot. To this 
end we write the particle current emerging from the left 
(right) reservoir in the form 



iV T ^ = ^(> 'ct, ,c 



:(£> 



L(R) ~ dV^ > Hp) k (p)> 
k( P ) 



f^(R)M • ( A1 ) 



Likewise, the energy current emerging from the left 
(right) reservoir can be shown to be given by 



Er.r-a^ - ^7(2^£fc(p) c fef c ) c fc(p)) 



du> 



; w2: l(r)( w ) 



L(R) - dt \Z_^ c Hp) L -k(p) u k(p)/ - / 27r ^^-L(R)^ 

Hp) 



The Green function calculation yields 



(A2) 



X L(R)( W ) = -* r L(,R)M(G o(u;) 

-/L(R)M[G a oH-GSoH]) , (A3) 

where the superscripts <, a, and r, denote the lesser, ad- 
vanced, and retarded Green function, respectively. The 
Fermi distributions, /l(r), are given in Eq. (|15[) . and the 
partial widths of the resonance level, r L ( R ), in Eq. (|16p . 
The Green function on the dot is calculated up to sec- 
ond order in the electron-phonon coupling 7^ One finds 






(w) = Goo M (^p M + Sf H) Gg ( W ) , ( A4) 
G00M = (« - e - E[(w) - <5e P - SpM)^ . (A5) 



where 



Here, 



<5e P = 2z 



•7 2 /" da V<f ^ 



(A6) 



is the polaron energy shift, where Qq Q (uj) is the lesser 
Green function on the dot in the absence of the coupling 
with the oscillator, 

G ™ {oj) ~ l (« - e Q y + (r( w )/2)» ' (A7) 

with T = Tl + r R . Here we have ignored a possible 
shift in the resonance energy due to the coupling with 
the leads, since it is not expected to play a significant 
role. 

As is seen from Eqs. (|A4[) and (|A5[) . the self energy on 
the dot includes two contributions. The first, £j, is due 
to the coupling with the leads, 



(A8) 



The second contribution to the self-energy results from 
the interaction with the phonons, and in second-order in 
7 reads 

yr (, A -,V f ^ ( (1 + WoM - NQ<W) 



and 



y< 



sr(w) = --^r L ( w ) + r R ( w )J , 

E<M =*(r L ( W )/ L (w) + r R (w)/ R (a;) 



27r V w — w — w' + i0+ 
A^>(q/)-(l + A)g<(t/h ^ 
cj + cj — ui' + i0 + / 



(w) = 7 2 (^^(w - w ) + (1 + A r )a < (w + w )) , 

(A10) 



where A denotes the phonon population. The lesser 
Green function Q < is given in Eq. (|A7|) . and the greater 
one, Q > , is given by the same expression with the distri- 
butions /l.r replaced by /l,r — 1. 

Inserting the expressions for the Green function Goo 
into Eq. (|A3[) . one finds that 2j/ R ) (w) can be written as a 
sum of two terms, one arising from the elastic transitions 
of the transport electrons, and the other coming from the 
inelastic ones, 

^L(R) (W) = 2L(R) M + 4Tr) ( W ) ■ ( A11 ) 

The elastic-process contribution is 

Z£(w) = |G5 (w)| 2 r L (w)r R ( w )[/ R ( w ) - f h (u)] , (A12) 

while the inelastic one is proportional to the strength of 
the electron-phonon coupling, 



2TV) -7 2 r L Ml^ (^)| 2 |ao r o(^-^o)l 2 

x 53 r a («-« )[JV/ a (w-w )(l-/ L (u;))-(l + iV)/ I( (w)(l-/ a («-Wo))] 



a=L,R 
,,2-p /-, .M/->r /, .\\2\nr 



- 7 ^r L H|^ HHao r o(^ + ^ )l 2 

x ^ r a ( W +« )[JV/ L (w)(l-/ a («+w ))-(l + iV)/ a (« + w )(l-/ L (a;))] ) 



a=L,R 



r 



(A13) 



where 5g , the retarded Green function in the absence 
of the coupling to the vibrational modes, is given in 
Eq. (|T9| . Since 2r is obtained from Eqs. (|A12[) and 
(IA13|) upon interchanging L with R, it is easy to see 
that the elastic-process parts of both the particle and 
the energy currents arc conserved (this is so because 
lf;{io) +Zfl(o;) = 0). The consideration of the inelastic- 
process part is a bit more delicate. By changing the in- 
tegration variables one finds that 



=^irV) = 7 2 / ^isjoKiPra"-)! 2 



x(r L ( w+ )r L ( w _)( w ;-wi)F LL (u;) 

+ w+r L (w + )r R (u;_)F LR (u;) 

-^r R K)r L ( w _)^ RL H) , » = oori, (Ai4) 



where F aa i is given in Eq. (|21[) . Here, w± = w ± 
luq/2. Hence, the inelastic-process parts of the parti- 
cle current (for which s = 0) are also conserved, i.e., 

V)] = 0. Using Eqs. (|A"I|l . 
produces Eq. (fTT|) for the 



f(du/2ir)[lt el (u) 



X R 



(|AT2|) . and (|A"T4)) in Eq. 
charge current. 

On the other hand, the energy current carried by the 
electrons alone is not conserved, since [using Eq. (|A14[) 
with s = 1] 

^[2L nol M +4 ncl M] = -o7 2 / ^l%K)l 2 
xieoVcOl 2 E r Q ( w+ )r Q ,( w _)F QQ ,( w ) . (Ai5) 

a,a'— L,R 



This result, in conjunction with Eqs. (|10|) and (|A2|) . 
leads to Eq. (|2U)) . Finally, the net energy current carried 
by the electrons [see Eq. ©] is obtained by using Eq. 
(|A14[) in Eq. (|A2[) (and the corresponding equation for 
E R ). This yields Eq. f[22]). 

In the case of the "floating" molecule, which is not 
coupled to a heat bath of its own, it is straightforward 
to show [using the Hamiltonian (fTT|) ] that the rate of 
change of the phonon population^ is given by minus the 
right-hand side of Eq. ()A15|) divided by ujq; i.e., 



27 



u>[r^\ u ) 



T-incl 
Z R 



M] 



dN 
] ~dt 







(A16) 



yielding Eq. Q for the energy conservation, with 
E-p = uj dN/dt. However, the phonon population of 
a floating molecule will arrange itself according to the 
chemical potentials and the temperatures of the elec- 
tronic reservoirs^ Consequently at steady-state dN/dt 
will vanish, implying that El + £r = 0, such that the 
energy current of the electrons is conserved. 

Finally we estimate the rate of decay of the vibration 
population due to the coupling with the electrons in the 
leads, when the latter are at thermal equilibrium. In 
diagrammatic language that rate is given by dressing the 
"phonon" line with an electron bubble, i.e. 



dN 
~dt 



duT(uj + ) Ma r(w_) 2 

yoo^+n — o — |yoo( w -) 



2tt 2 ^ uuv +/l 2 

k (jV/(w_)[l - /K)] - (1 + N)f(u + )[1 - /(w_)]) . 

(A17) 

At zero temperature the last factor in the integrand limits 
it to the range \u — fx\ < cj /2, leading to the rate 



(a,o/[27r])( 7 /r) 2 



(A18) 



when r 3> loq, as mentioned in Sec. HI In general, the rate 
is a non-monotonic function of the ratio tuo/T, reaching a 
maximal value when these two energies are comparable. 
The result of Eq. (|A18|) can be qualitatively obtained in 
a more elementary, but equivalent, fashion using third- 
order perturbation theory. Consider the decay of the 
excited vibrational mode with a T = Fermi gas on 
each lead (having a density of states vq). For simplic- 
ity we take the resonant case and assume r 3> ujq to 
start with. For the decay with the left lead, the first in- 
termediate state has an electron from that lead go into 
the molecule with an amplitude V£ and an effective en- 
ergy denominator Tl/2 (due to being on resonance), in 
the second intermediate state, the vibration is deeexcited 
(amplitude 7) and the electrons stay put. For r ^> luq, 
the energy denominator is again approximately Tl/2. In 
the final state, the electron goes back to the same lead, 
with an energy loq higher that the one it started from, 
and with an amplitude Vl- The total amplitude for this 
process is V^7Vl/(Fl/2) 2 . Finally, taking the absolute 
square of the amplitude and multiplying by 27Wo, we get 
the goldcn-rulc rate for this decay. Multiplying by the 
number of such processes, vqu>o, and summing over the 



(equivalent) leads, we get Eq. (|A18|) in order of magni- 
tude. Clearly, in the opposite case r <C u>q, one T in 
the denominator is replaced by ojq. A surprising feature 
of this result, which must be pointed out, is its decrease 
with r, the rate to get into/from the molecule from/into 
the leads. This is quite counterintuitive, but it is what 
the quantum-mechanical calculation tells us! Formally, 
this is due to Green function at the resonance having iT/2 
as its denominator, meaning physically that the width of 
the resonance sets the scale for the "closest approach" to 
it. 
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